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ABSTRACT 



Characteristics of sea and swell incident at a permeable rubble- 
mound breakwater located in Monterey Harbor, California, are resolved 
into reflected and transmitted components. The wave characteristics 
are studied by analyzing synchronized wave records of three underwater 
sensors judiciously placed, two to seaward and one to landward of the 
breakwater. Power spectra and cross spectra are calculated for various 
characteristic sea states selected from three months of observations. 
Amplitude and phase are determined for the spectral wave components 
comprising the partial standing wave phenomena. This study was unique 
in that it entails experiments conducted in the field on a prototype 
structure in the natural environment. Transmission has been studied 
heretofore exclusively with scale models. The coefficients of trans- 
mission are considerably less for the prototype than predicted by a 
Corps of Engineers model study of the Monterey Breakwater. The differ- 
ences are apparently related to more wave energy scattering and dissi- 
pation due to turbulence than predicted by the model. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION 10 

A. BACKGROUND 10 

1. Monterey Harbor 10 

2. Breakwater Characteristics 13 

B. PROBLEM AND APPROACH 16 

C. OBJECTIVE 18 

II. THEORETICAL CONSIDERATIONS 20 

A. WAVE DESCRIPTION 20 

B. WAVE AMPLITUDES 20 

C. WAVE REFLECTION AND TRANSMISSION COEFFICIENTS 25 

D. WAVE ENERGY DISSIPATION 25 

III. ANALYTICAL PROCEDURE 27 

A. SEAWARD WAVE DATA 27 

B. HARBOR WAVE DATA 27 

C. COMPUTATION OF POWER SPECTRUM 30 

D. COMPUTATION OF CROSS SPECTRUM 31 

IV. INSTRUMENTATION 36 

A. SENSOR LOCATION 36 

B. SEAWARD SENSORS 38 

C. HARBOR SENSOR 39 

D. WAVE RECORDERS 39 

V. ANALYSIS OF RESULTS 42 

A, INTERPRETATION OF POWER SPECTRA 42 

1. Measured Power Spectra 4 ^ 

2 . Incident and Reflected Sea Surface Power Spectra . . 58 



3 



B. INTERPRETATION OF CROSS SPECTRA 59 

C. TRANSMISSION, REFLECTION, AND ENERGY 

DISSIPATION COEFFICIENTS 62 

D. COMPARISON WITH MODEL STUDY 73 

E. ERROR ANALYSIS 77 

1. Data Reduction 77 

2. Instrument Location 78 

3. Numerical Analysis 78 

VI. CONCLUSIONS AND SUMMARY 80 

REFERENCES 83 

INITIAL DISTRIBUTION LIST 85 

FORM DD 1473 87 



4 



LIST OF FIGURES 



Figure Page 

1 Monterey Harbor H 

2 Monterey Bay 12 

3 Monterey Harbor Refraction Diagrams 14 

4 Monterey Breakwater 15 

5 Monterey Breakwater Structure 17 

6 Transmission through a Permeable Breakwater 21 

7 Geometric Relationships for Conversion from 

Curvilinear Coordinates (R,9) to Rectilinear 
Coordinates (Y,t) 29 

8 Linear Interpolation 29 

9 Seaward Sensor Positions 37 

10 Harbor Sensor Position 37 

11 Wave Sensor Mounts Before Study (above) 

and After Study (below) 40 

12 Measured Power Spectra, High Tide, 

17 December 1970 (1330) 43 

13 Phase Spectrum, High Tide, 17 December 1970 

(1330) 43 

14 Incident, Reflected, and Transmitted Power 

Spectra, High Tide, 17 December 1970 (1330) 44 

15 Coherence Spectrum, High Tide, 17 December 

1970 (1330) 44 

16 Measured Power Spectra, High Tide, 

1 December 1970, (0800) 45 

17 Phase Spectrum, High Tide, 1 December 1970 

(0800) 45 

18 Incident, Reflected, and Transmitted Power 

Spectra, High Tide, 1 December 1970 (0800) 46 

19 Coherence Spectrum, High Tide, 1 December 

1970 (0800) 46 



5 



20 Measured Power Spectra, Low Tide, 

1 December 1970 (2000) 47 

21 Phase Spectrum, Low Tide, 1 December 1970 

(2000) 47 

22 Incident, Reflected, and Transmitted 
Power Spectra, Low Tide, 1 December 1970 

(2000) 48 

23 Coherence Spectrum, Low Tide, 

1 December 1970 (2000) 48 

24 Measured Rower Spectra, Low Tide, 

1 December 1970 (2330) 49 

25 Phase Spectrum, Low Tide, 1 December 1970 

(2330) 49 

26 Incident, Reflected, and Transmitted 
Power Spectra, Low Tide, 1 December 1970 

(2330) 50 

27 Coherence Spectrum, Low Tide, 

1 December 1970 (2330) 50 

28 Measured Power Spectra, Low Tide, 

2 December 1970 (2030) 51 

29 Phase Spectrum, Low Tide, 2 December 

(2030) 51 

30 Incident, Reflected, and Transmitted 
Power Spectra, Low Tide, 2 December 1970 

(2030) 52 

31 Coherence Spectrum, Low Tide, 

2 December 1970 (2030) . . . 52 

32 Measured Power Spectra, Low Tide, 

14 December 19 70 (1600) 53 

33 Phase Spectrum, Low Tide, 14 December 1970 

(1600) 53 

34 Incident, Reflected, and Transmitted 
Power Spectra, Low Tide, 14 December 1970 

(1600) 54 

35 Coherence Spectrum, Low Tide, 

14 December 1970 (1600) 54 

36 Analog Trace Illustrating Partial 

Standing Wave Pehnomenon 57 



6 



37 Coefficients of Reflection and Transmission, 

High Tide Conditions ^4 

38 Coefficients of Energy Dissipation, High 

Tide Conditions 65 

39 Breakwater at High Tide, 1 December 1970 66 

40 Breakwater at High Tide, 17 December 1970 67 

41 Coefficients of Reflection and Transmission, 

Low Tide Conditions 69 

42 Coefficients of Energy Dissipation, 

Low Tide Conditions 70 

43 Breakwater at Low Tide, 14 December 19 70 71 

44 Transmitted Wave Heights for Prototype 

Study and Hydraulic Model Investigation 76 



7 



LIST OF TABLES 



Table 

I 

II 



Monterey Breakwater Wave Data from Field Study 

« 

Monterey Harbor Wave Transmission Data from 
Hydraulic Model Investigation 



Page 

63 

• • 

74 



8 



ACKNOWLEDGEMENTS 



The author wishes to express his sincere appreciation to Dr. Edward 
B. Thornton of the Department of Oceanography, Naval Postgraduate School, 
under whose direction this thesis was written. His encouragement, 
criticism, recommendations, and experience were invaluable in the 
completion of this research. Special thanks are also given to Professor 
Warren C. Thompson for critically reading the manuscript. 

In addition, the author is indebted to Mr. Robert C. Smith, Research 
Technician with the Department of Research Administration, for his 
valuable assistance in calibrating and preparing the wave sensors for 
operation, and to Mr. Jack C. Mellor, Oceanographer with the Department 
of Oceanography, for his technical advice and assistance in installing 
the wave sensors. 



9 



I. INTRODUCTION 



A. BACKGROUND 

1. Monterey Harbor 

Monterey Harbor (Figure 1) is located at the southern end of 
Monterey Bay about 75 nautical miles south of San Francisco, California. 
The Bay is a large symmetric, cres cent -shaped , open embayment approxi- 
mately 22 nautical miles long from north to south, nine nautical miles 
wide, and is slightly closed at its extreme northern and southern ends 
by two stubby land masses (Figure 2) . A unique feature of the bay is 
the deep Monterey Submarine Canyon which divides the bay roughly into 
two halves. The submarine canyon is one of the largest in the world, 
and it exerts significant influences on both directional and height 

characteristics of waves within the bay. 

Until the 1940's the harbor's primary function was to provide 

an anchorage haven for fishing fleets. However, due to the reduced 
fish catch in recent years, the activities of the harbor have been 
largely reoriented to recreational boating and tourist attraction fea- 
tures including sports fishing. In 1934, a rubble-mound breakwater was 
constructed to provide harbor protection for the expanding fishing 
industry. It soon became evident that this structure did not provide 
sufficient mooring area for the growing number of fishing and small 
craft, nor did it provide adequate protection from wave and surge action. 
Current plans include construction of additional breakwater structures 
to provide a larger and safer anchorage for small craft. 

Long-period waves occurring in the bay at times cause a sub- 
stantial surge in the harbor. In addition, the harbor is exposed to 
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(View looking west) 




(View looking northwest) 



FIGURE 1 
Monterey Harbor 
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Monterey Bay 
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sea and swell arriving from outside the bay from the deepwater directions 
between west and north-northwest, and on rare occasions to local winds 
from the north within the bay [Bixby, 1962]. Refraction diagrams 
(Figure 3) obtained from the San Francisco District Corps of Engineers 
show that waves entering the bay from the deepwater directions from 
approximately west to north-northwest (260 deg to 332 deg) are almost 
unidirectional in the vicinity of the harbor. These waves are often of 
sufficient size to damage harbor facilities as well as fishing boats 
and pleasure craft. There have been several studies concerned with 
wave and surge effects on the bay and the harbor which have led to 
analytical solutions and hydraulic model investigations [Hudson, 1949; Wilson, 
Hendrickson and Kilmer, 1965; Chatham, 1968; Bucci and Whalin, 1969]. 

2. Breakwater Characteristics 

The Monterey structure is a typical permeable, rubble-mound type 
preferred for coastal harbors because it is economical, adaptable to any 
reasonable water depth, suitable on nearly all foundations, and readily 
repaired. Rubble-mound breakwaters are more or less a heterogeneous 
assemblage of natural stones of different sizes and shapes either dumped 
at random or placed in courses. Side slopes and armor stone sizes are 
designed to effectively resist the expected wave action. The rubble-mound 
type is used extensively throughout the United States and almost exclu- 
sively on the Pacific Coast. 

The first Monterey Harbor breakwater construction was a 1300-foot 
rubble-mound structure completed in early 1934. A 400-foot extension 
was completed late the same year (Figure 4) . The crest is at elevation 
plus ten feet (MLLW) and is 15 feet in width. The harbor side slope is 
one and one-half horizontal on one vertical for the full length of the 
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'lonterey Harbor Refraction Diagrams 
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FIGURE 4 

Monterey Breakwater 
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breakwater. The seaward side slope is one and one-half horizontal on 
one vertical for the original 1300 feet. The 400-foot extension has 
a seaward side slope of one and one-half on one vertical from crest to 
elevation minus 16 feet (MLLW) , and one and one-quarter horizontal on 
one vertical below that depth. The structural adequacy of the break- 
water has been substantiated by the fact that no replacement of stone 
or other preventative maintenance has been required since it was built, 
and it has remained a stable configuration. 

The Monterey breakwater consists of an interior section, or core, 
of assorted sizes of stone, gravel, and fines, protected by a ten-foot 
layer of coarse, angular-shaped armor cover components (Figure 5). The 
exterior A stones, or primary cover layer, placed to attain an inter- 
locking fit, are approximately ten tons per rubble component for the 
original 1300 feet, and four to twelve tons each for the 400-foot 
extension. The original 1300 feet has a cement cap roadway laid in 
conjunction with the construction of a wharf and mooring dolphins on 
the harbor side. The breakwater has never been grouted. 

B. PROBLEM AND APPROACH 

When sea and swell impinge upon sloping permeable rubble-mound 
structures, they may break, and thereby dissipate energy by turbulence, 
they may be reflected to produce a partial standing wave by inter- 
ference with the incident waves, and they may be transmitted through 
or partially over the structure with energy loss due to turbulence 
within the breakwater. The resulting interplay of forces developed by 
the wave-induced water motions and the action of the rubble units is 
extremely complex due to the variable reflective and friction properties 
of the sloping face and the general characteristic roughness of the 
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FIGURE 5 

Monterey Breakwater Structure 
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interior of the permeable structure. Attempts to determine by 
theoretical analysis the hydrodynamic characteristics related to 
stability and transmission have not been productive. 

The problem of determining wave transmission was surmounted in 
this study by application of the principal of energy conservation to 
statistical wave observations. This was done by measuring the wave 
motion in the field with wave sensors judiciously placed on both sides 
of a permeable breakwater. Wave components were computed through 
spectral analysis of the individual wave records. Reflected and 
incident wave components were determined through application of linear 
wave theory for each frequency bandwidth. The transmitted wave component 
was computed from measurement of the progressive wave in the lee of the 
breakwater. The energy dissipated was then quantitatively deduced by 
determining the difference between the incident wave energy and the 
reflected and transmitted wave energy. 

Sample wave records from this particular wave field were considered 
accurate estimates of an ergodic, stationary, random process. Con- 
version of these finite sample time histories to a spectral density 
(frequency) representation was accomplished through computer analysis. 
Fourier analysis was used to break down the random sea process into 
simple sinusoidal harmonic components whose behavior was analyzed by 
the use of classical linear theories of wave motion. 

C. OBJECTIVE 

The objective of this investigation is threefold, and can be con- 
veniently separated to include the following areas of interest: 

1. Investigate the feasibility of determining the incident and 
reflected components on the incident side of a permeable breakwater 
through selective placement of wave sensors in the field. 
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2. Consider the wave -transmission problem as a stochastic process as 
opposed to consideration of monochromatic waves; utilize spectral analysis 
to resolve the random sea environment into a continuum of wave components 
described by linear wave theory for a specified range of frequencies; 

and compute transmitted and reflected wave coefficients and the energy 
dissipation of the breakwater from equations derived from linear wave 
theory, 

3, Compare the results of this field study with hydraulic model 
transmission tests of the combination tribar, rubble-mound breakwater 
proposed for Monterey Harbor by Davidson [1969]. 
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II. THEORETICAL CONSIDERATIONS 



A. WAVE DESCRIPTION 

The sea and swell selected for this study arrived incident at 
Monterey Harbor approximately normal to the breakwater (Figure 3). The 
areas adjacent to the breakwater are fairly uniform in depth. It is 
assumed that these random sea and swell waves can be analyzed as the 
result of the superposition of a number of component sinusoidal (Airy) 
wave trains with their own characteristic height, frequency, direction, 
and random phase. For incident waves of low steepness at the depths 
involved in this study. Airy theory was shown by Dean [1968] to give a 
reasonably valid approximation. These considerations are extended to 
random sea theory. The linear wave theory problem for a partially 
reflected wave is presented below. The linear solutions are then super- 
imposed to obtain the spectrum of the actual sea surface. 

B. WAVE AMPLITUDES 

The incident and reflected wave amplitudes are computed by consider- 
ing two progressive waves moving in opposite directions comprising a 
system of imperfect or partial reflecting waves (Figure 6). Basic 
assumptions include: monochromatic and long-crested waves, wave 

direction normal to the breakwater, and a horizontal seabed. Assuming 
linear theory, the incoming and reflected components of the wave surface 
elevation, n, can be superimposed as follows, 

= n-]-+n R = a^cos (kx^t+dj) + a R cos (kx^at+6^) (1.1) 

n 2 = = SjCos (kx 2 ~ot+6 I ) + a R cos (kx 2 +at+6 R ) (1»2) 
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WAVE DIRECTION 




FIGURE 6 

Transmission through, a Permeable Breakwater 
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where , 



H = height of wave, ft. 

u 

a = amplitude of wave (-j ) , ft. 

L = wave length, ft. 

k = wave number (^-) , ft ^ . 

Li 

T = wave period, sec. 

2tt ”1 

a = angular frequency (— ) , sec 
6 = phase angle, degrees. 

Subscripts 1, 2, and 3 denote quantities corresponding to wave sensors 
numbered 1, 2, and 3; and I, R, and T refer to incident, reflected, and 
transmitted wave components. 

The measured statistical values of sea surface elevation from the 
wave sensors are , 
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cos 


(a t 


- V 


(1.3) 


= A 2 


cos 
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- V 


(1.4) 


= a T 


cos 


(at 


- S 3> 


(1.5) 



where , A^ y and a^ are statistical amplitudes measured at wave sensors 
numbered 1, 2, and 3, respectively. 

Using trigonometric identities, equations (1.1) and (1.2) are 
expanded as , 



tu = [a_ cos (kx.. + 6_) + a_ cos (kx- + 6_)] cos at 
i I 1 I R 1 R (1.6) 

+ [a T sin (kx- 4- 6 ) - a_ sin (kx- + 5^)] sin at 
1 ± 1 K X K 

n 0 = [a T cos (kx_ + 5 ) + a D cos (kx~ + 5 )] cos at 

Z L Z L K Z K. 7) 

+ [a^ sin (kx£ + 5^.) - a^ sin (kx^ + 6^) ] sin at 



Similarly, trigonometric expansion of equations (1.3) and (1.4) 
yield , 
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( 1 . 8 ) 



(cos ot cos 6^ + sin at sin 6^) 

n 2 = A 2 ^ COS at C0S 5 2 + sin at sin ^2^ (1.9) 

Equating coefficients of cos at and sin at of the measured and 
theoretical values of and introduces the relationships. 
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Rewriting equations (1.10) and (1.11) using trigonometric identities 
yields , 
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Rearranging terms , 

cos 6^ = cos kx^ (a^ cos 6^. + 

- sin kx^ (a^. sin 6^. + 

A^ sin 6^ = cos kx^ (a^. sin 6^. - 

+ sin kx^ (a^ cos 6^. - 

Similarly rewriting equations (1.12) 
expansion yields , 

A 2 cos 6 2 = cos kx 2 (a^. cos + 

- sin kx 2 (a^. sin 6^ + 
A 2 sin 6 2 = cos kx 2 (a^. sin 6^. - 

+ sin kx 2 (a^ cos 6^ - 



^ cos 6 r ) (1.16) 
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Judiciously choosing the origin of x at wave sensor No. 1 sets 
quantity x^ equal to zero. Further, since phase angles are relative 
quantities, 6^ is set equal to zero which redefines 6^ as the phase 
difference, henceforth denoted as 6. Thus equations (1.16) and (1.17) 
are further simplified to yield, 

A^ = a^ cos 6^ + a R cos 6 R (1.20) 

0 = a^ sin 6^ - a R sin S R (1.21) 

Substitution of equations (1.20) and (1.21) into (1.18) and (1.19) 
yields , 
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( 1 . 22 ) 
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Solving for sin 6 and cos 6 in equations (1.22) and (1.23) and 
K K 

dividing yields , 



S R = tan 



-1 



cos kx^ - A^ cos 6 
A^ sin kx^ - A 2 sin 6 



(1.24) 



Solving for a in equation (1.22) in terms of 6 gives, 
R R 



-A 2 cos 6 + cos kx 2 
^ — . 

R 2 sin kx„ sin 6, 



(1.25) 



2 ““ ~R 

Utilizing equations (1.20) and (1.21) and dividing gives the solution 



for 6 in terms of 6 and a as, 
i R R 



5^ = tan 



-1 



a R Sl " { R 

A r a R cos 5 r 



(1.26) 



Solving for a^. from equation (1.21) gives the solution for a^ in 

terms of 6 , 6 , and a yields, 

1 K K 



a i = 






sin 6 



R 



sin 6. 



(1.27) 
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C. WAVE REFLECTION AND TRANSMISSION COEFFICIENTS 



The quantities and a^ have been evaluated in terms of known or 
measurable quantities A^, k^, k, and 5. The reflection coefficient, 
and transmission coefficient, K^, are defined as, 



K = — 
R a„ 



(1.28) 



ar 

K = — 
T a x 



(1.29) 



and can therefore be evaluated using computed values a„, a T , and 

R X 



measured quantity a^,. 



D. WAVE ENERGY DISSIPATION 

The average power per unit length of wave crest can be expressed 
analytically as the average time rate of change of energy within an 
element of fluid through which the wave passes. The power dissipated 
by the breakwater can be computed by writing an expression equating 
the power incident on the breakwater, the power transmitted through 
the breakwater, and the power reflected from the breakwater. The 
equation for a power balance for both the seaward and the harbor side 
of the breakwater may be written as, 



a y a T a R 

Y "2 C G r ~ Y ~2 C G t ~ Y “T C G r = ? D 



(1.30) 



where , 

3 

Y = specific weight of sea water, 64.0 lb/ft . 

C = wave group velocity, ft /sec. 

Pp = power dissipated, ft-lb/sec per ft of wave crest. 

Assuming constant depth for simplification and rewriting the terms 
in equation (1.30) yields, 



25 



(1.31) 
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Using previously defined notation yields, 






(1,33) 



where , 



= coefficient of energy dissipation expressed as the fraction 
of the incident wave power dissipated. 



The simplifying assumption of constant depth is justified to seaward 
of the breakwater. The depth changes from 47 feet at the seaward 
sensor (No. 1) to 42 feet at the toe of the breakwater (Figure 4). The 
depth decreases from an average of 46 feet for the two seaward sensors 
(No. 1 and No. 2) to 39 feet at the harbor sensor (No. 3). The assumption 
of a constant depth bottom of 46 feet results in a maximum error for the 
transmitted energy of less than eight per cent for a 20-second wave, the 
error decreasing with decreasing period, and so will be used in the 
analysis. 
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III. ANALYTICAL PROCEDURE 



A. SEAWARD WAVE DATA 

The wave records to seaward of the breakwater are rectilinear 
analog traces measured with two identical sensors. These records were 
digitized at 100 points to the inch using a Calma Company Model 480 - 
mechanical digitizer and interpreted by program CONVERT [Lynch, 1970] 
using the CDC 6500 computer at Fleet Numerical Weather Central, Monterey. 
A record length of 2048 points was digitized with a sampling interval 
of 0.63 seconds for a time record of 0.358 hours or approximately 21.5 
minutes. This corresponds to approximately 250 oscillations of the 
average period wave which is generally considered a reliable base 
sample for analysis of ocean wave records. The sensors eliminate 
tidal, seiche, and other long-period influences through high-pass 
filtering by hydraulic means inherent in the pressure differential 
transducers (time constants 47 and 50 seconds). Possible electronically 
generated direct-current drift was removed by computing the linear trend 
and correcting the individual raw data points. 

B. HARBOR WAVE DATA 

The wave records from the harbor side of the breakwater are curvi- 
linear analog traces measured with a single sensor. Although a differ- 
ential pressure transducer of different design than the seaward wave 
sensors was used, high-pass filtering was similarly accomplished by 
hydraulic filtering (time constant 60 seconds). The raw data were 
corrected for linear trends using the same procedure as for the seaward 
data. 
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As previously stated, the harbor wave data were recorded on curvi- 
linear coordinates. The Calma 480 digitizer accepts rectilinear 
coordinates which resulted in additional problems in data reduction. 
These records were analyzed at 100 points to the inch, but at a 
sampling interval of 0.80 seconds. A record length of 2048 points was 
digitized for a time record of 0.455 hours, or approximately 27.3 
minutes. The unadjusted amplitudes and time increments obtained from 
the digitizer and computer program CONVERT were recorded in incremental 
paths every 0.01 inch of stylus travel. The rectilinear coordinates 
from the digitizer were then converted to the approximately correct 
curvilinear coordinates (t, R9) by employing geometric relationships 
as follows (Figure 7) , 



+ R 




and 



9 



i 



-1 

tan 




(3.1) 



where , 



= time in curvilinear coordinates. 

R = radius of curvature for arcs of constant t^. 

y^ = the ordinate on the cartesian scale. 

x_£ = time scale in equally spaced increments of t on the 
rectangular cartesian scale. 

The adjusted heights are then computed as, 

y ± = R9. (3.2) 

where is now the adjusted wave surface elevation value. Since the 
space increment is 0.01 inch, the adjusted time, t^, is, 

t ± = 80 x. (3.3) 
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Ua 




FIGURE 7 

Geometric Relationships for Conversion from Curvilinear Coordinates 
(R,9) to Rectilinear Coordinates (y,t). 




Linear Interpolation 
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A linear interpolation scheme is then performed to match the seaward 



data time increments. The computer checks each successive time value 
of the harbor wave record and interpolates to the appropriate increment 
of the uniform sampling interval. The slope of the line (Figure 8) is, 



Aj, frl - ?!-!> 

“ ' <‘i - C i-1> 



(3.4) 



The horizontal spatial scale for the interpolated values is, 

x i = (i-l)A (3.5) 

where A = the uniform sampling interval. 

This final interpolation to a uniform sampling interval, A, is a 
critical step in curvilinear to rectilinear coordinates conversion. 
Although not a requirement in the analysis , all the wave data are then 
of the same sampling interval of 0.63 seconds of wave time. 



C. COMPUTATION OF POWER SPECTRUM 

The analysis of the digital records is made in accordance with the 
method of Blackman and Tukey [1958] with the computations performed on 
an IBM 360 computer at the Naval Postgraduate School. This particular 
method of spectral analysis involves the determination of the auto- 
covariance function, which is first smoothed, and then the power 
spectrum is computed through Fourier transformation. The auto- 
covariance function of the wave record time series, X^(t) is expressed 

as the average lagged product, and is defined as, 

T/2 

$ u (t) = ^ t / X l (t) V t+T > dt (3.6) 

-T/2 

where t = time lag. 

Due to the finite nature of the given continuous record, the auto- 
covariance function is expressed within set limits for all values of t. 
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This finite expression is called the apparent auto-covariance function 
and is written as , 

(T n - I x I )/2 

♦„(t) - T 1 | T | f 3^(0 X 1 (t+x) dt (3.7) 

n -d n - |x|)/2 

where T = the length of the record, 
n 

A filtering function, called a lag window, is introduced to restore 
the features of the true spectrum. A Parzen lag window is used here 
and is defined as, 

D (t) = 1 - (% (1 - |), k = (3.8) 

P m 

where m = the total number of lags. 

The Parzen window function is always positive, i.e., it has no 
negative side lobes. This window is particularly useful in calculating 
cross-spectra as it tends to eliminate any numerical instability prob- 
lems that may arise in the analysis procedure. The modified, apparent 
auto-covariance function, <t>p(x), is then given as, 

4> p ( T > = D p (x) ^(x) (3.9) 

The smoothed spectral estimates are determined by taking the cosine 

transform of the modified, apparent auto-covariance function. This 

final representation is called the power spectrum and is expressed as, 

00 

♦ (a) = f <}> (t ) cos axdx (3.10) 

11 2n J P 

—00 

D. COMPUTATION OF CROSS SPECTRUM 

For the problem at hand the cross spectra of pairs of stationary, 
random times series contain two different but related and extremely 
useful forms of information: the phase spectrum and the coherence 
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spectrum. The phase spectrum measures the phase difference between the 
two processes, and the coherence spectrum is a measure of the correlation 
between the two processes, at each spectral frequency band. 

By definition, the cross spectrum, $ 22 * t * ie Courier transform of 
the cross-covariance function with both odd and even functions of their 
respective arguments, and can be written as, 

00 

$l 2 ( a ) = Yn j 4*12 ( T )( cos ox - i sin at) dt (3.11) 

— oo 



where = the wave record time-series for sensor No. 1. 

= the wave record time-series for sensor No. 2. 

<{>^2 = the cross-covariance function between X^ and X^ given by. 



T/2 

♦ 12 (t) = f X 1 (t) X 2 (t+r) dt (3.12) 

-T/2 



In the actual computer analysis, the cross-covariance function is 
modified by the Parzen lag window, but for simplicity in terminology, 
the apparent, modified cross-covariance function will be referred to as 
the cross-covariance function. 

The cross spectrum as a complex-valued function of a can be 
written as 



* 12 (t) 



■Hi 






cos axdx 



00 



4>^ 2 (x ) sin axdx 



(3.13) 



This spectral function combines both the cosine and sine transforms of 
2 (x ) anc ^ can resolved into odd and even complements. Since the 
cosine (sine) transforms of an odd (even) function is zero, the cross 
spectrum can then be expressed in this simplified form. 
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*12 (o) = h ) f *even (t) COS ° xdx-i / *odd (x) Sin aX£ 4 (3 ’ 14) 



A useful property of the cross-covariance function which will be used 
extensively later in the analysis can be shown by replacing x by -x , and 
then making the substitution 6 = t-x to obtain. 



, -v _ lim 

4*12 T-h*> 



T/2 • 

| ^(5+t) X 2 (S)dS 



(3.15) 



-T/2 

Since the variable of integration here is a dummy variable, a valid 
expression of the cross-covariance as an even function is, 

<f>l 2 (“ x ) = 4> 21 ( t ) (3.16) 

At this stage of the analysis, useful functions A.^(x) and B^ 2 (x) 
are defined by. 



Ai 2 ( x ) = 



B 12 (t) = 



<f> 12 (T) + *21^ 



<t>l 2 ( x ) ~ ^("O 



(3.17) 



Replacing x by -x in each expression and making use of the relationships 
of (3.16) yields the expressions, 

A 12 (-x) = A 12 (x), and 



(3.18) 






where A^ 2 (x) as an even function and B^ 2 (x) is an odd function. 
Since, by algebraic manipulation, <}>^ 2 (x) can be resolved as. 



4>1 2 ( x ) = A 12 (x) + B 12^ X ^ 



(3.19) 



the functions A^ 2 (x) and B 12^ X ^ are > respectively, the desired even and 
odd parts of the cross-covariance function. Expression (3.14) can now 
be written as, 
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(3.20) 



$ 12 (a) = j J A 12 ( T ) cos 0TdT-i J B 12 ( T ) sin aTdx f 

l, — <x> —CO J 



The real part in the above expression defines the co-spectrum, 
and the negative of the imaginary part defines the quadrature spectrum, 
Ql 2 (o). The cross spectrum is then written in simplified complex form 
as , 

*12 (a) = C 12 (a) “ 1 Q 12 (a) (3.21) 

It should be noted at this stage of the analysis that the co-spectrum 
is the cosine transform of the even part of the cross-covariance function. 
Analogously, the cosine transform of the auto-spectrum is the power 
spectrum of the record. The power spectrum is defined as showing the 
contribution to the total variance of a single record in the indicated 
bands of frequency, and has the same significance as the even part of 
the cross-covariance function; the evenness expresses the cross-covariance 
of the in-phase (zero degree) components or reversed phase (+ 180 degrees) 
components between the two records. 

The quadrature spectrum, with a similar argument, can be shown to 
represent the amount of power due to components which are out of phase 
by a constant amount (+ 90 degrees). 

The spectral phase difference, is defined as the argument 

of the complex conjugate of the cross spectrum and is written as, 

-1 ^ 12 ^ 

e i2 <0) * (3 - 22) 

The spectral phase difference represents the average phase angle by 
which the cross-correlated components of lead those of X^ in each 
spectral frequency band. 



The coherence is extremely important as a meaningful estimator of 
the creditability of the phase spectra in each spectral frequency band. 
Not all the frequency components from the two synchronized random 
records may be cross-correlated , contrasted with the perfect correlation 
(at zero lag) of the components of the individual power spectra. This 
comparison can be exhibited by a ratio of the cross spectrum to the 
square root of the individual auto-spectra and is written as, 



R (a) = 




This coherency function can then be expressed quite simply as , 



R (a) 



c * 2 (a) + Q 2 (a) 
$ ll (a) $ 22 (a) 



1/2 



< 1 



This ratio plays the role of a correlation coefficient defined for each 



spectral frequency band. 
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IV . INSTRUMENTATION 



A, SENSOR LOCATION 

The wave data obtained were recorded from three wave sensors placed 
in a line normal to the breakwater (Figure 4) • The distance between the 
two seaward sensors was pre-planned at 135 feet, which is one-quarter 
wave length of a 14-second wave in 46 feet of water. This period 
coincides with the most frequently occurring wave period for the months 
of November to January, which is approximately 14 seconds as determined 
from sea and swell data for incident waves in the vicinity of the harbor 
[Wilson, Hendrickson, and Kilmer, 1965]. The actual seaward sensor 
separation was determined by transit and range lines to be 130 feet. 

The seaward sensors (No. 1 and No. 2) are at depths of 45 and 47 feet, 
for an average depth of 46 feet (MLLW) . The locations of the sensors 
with respect to the breakwater can be seen in Figures 9 and 10. 

Although not imperative in the analysis of the partial standing wave 
phenomenon, the separation between sensor No. 2 and the breakwater face 
at MLLW was established at 135 feet in an effort to observe the standing 
wave under ideal conditions. Sensor No. 2 was approximately 135 feet 
from the breakwater face at MLLW, and approximately 100 feet horizontally 
from the base of the breakwater. This separation is considered necessary 
for wave records to be uneffected by the turbulence of the wave upwash 
or other related effects at the face of the breakwater. The horizontal 
separation between sensor No. 3 (harbor) and the breakwater face at MLLW 
is approximately 80 feet which is necessary for a balance between minimum 
diffraction and negligible turbulence effects. Location of the harbor 
sensor is such that diffraction may be considered negligible. 
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BREAKWATER 










SENSOR #2 



SENSOR#1 



FIGURE 9 

Seaward Sensor Positions 




FIGURE 10 

Harbor Sensor Position 
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B. SEAWARD SENSORS 



The two sensors used to seaward of the breakwater are identical 
Marine Advisors Model A-2b instruments modified from the basic Frank 
E. Snodgrass Mark IX design. The A- 2b sensor consists of two oil-filled 
chambers connected through a capillary tube and, in a parallel flow 
path, a strain-gauge differential pressure transducer. One chamber is 
exposed to sea pressure with fluid separation accomplished by a 
neoprene diaphragm. The second chamber serves to introduce compliance 
into the system by means of a spring-loaded, air filled bellows which 
allows flow through the capillary to slowly equalize the pressure 
between the chambers. Rapid pressure fluctuations, such as those due 
to wind-driven waves and swell, result in an insignificant amount of 
flow through the capillary, and a pressure differential is registered 
between the chambers. Slow variations in pressure, such as those due 
to tides and seiches, cause sufficient flow during intervals of time 
much shorter than their period to prevent a significant pressure dif- 
ferential from being maintained across the transducer. When in this 
configuration, the device becomes a hydraulic analog of an electrical 
RC high-pass filter circuit. The two A- 2b sensors are hydraulically 
filtered at time constants of 47 and 50 seconds and produce a uniform 
flat response curve for sea and swell data. 

When placed on the sea floor, the A- 2b transducer senses the 
pressure signature of the waves present in the form of short-period 
deviations from the average pressure at the instrument. The inter- 
pretation of the pressure fluctuations is modified by a hydrodynamic 
attenuation factor that is a function of the depth at the instrument, 
elevation of the instrument above the bottom, and the period of the 
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waves being observed. The sensors are tripod mounted (Figure 11). 

The seabed in the vicinity of the seaward sensors is a relatively 
flat, clay surface with bottom contours generally straight and 
parallel to the breakwater. These bottom characteristics result in 
negligible shoaling and refraction in the vicinity of the seaward wave 
sensors . 

C. HARBOR SENSOR 

The single sensor used on the harbor side of the breakwater is a 
Frank E. Snodgrass Mark IX design. The sensor f s principle component 
is a Bourns differential pressure potentiometer which is used as the 
transducer. Hydraulic filtering at a time constant of 60 seconds, 
comparable in theory and sensitivity to the Marine Advisors A-2b 
sensors, is achieved in the Mark IX model. The harbor sensor is in an 
approximate depth of 39 feet at mean lower low water. Linearity and 
resolution of + 1 per cent of full-scale deflection is obtained through 
sensitivity control adjustment of the instrument. This provides a uni- 
form flat response curve for sea and swell data. 

D. WAVE RECORDERS 

The seaward wave data obtained using the Marine Advisors A- 2b 
sensors were recorded by a solid-state, modular-construction, Clevite 
Brush recorder. Wave data from both sensors were recorded simultaneously 
with two single-channel, plug-in modules on pressure-sensitive paper 
with a rectilinear trace presentation. This multi-trace design insures 
the accurate synchronization between both sensors required for cross 
spectra calculations. The chart speed used was 0.4 millimeter per 
second. 
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FIGURE 11 

Wave Sensor Mounts 

Before Study (above) and After Study (below) 
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Manufacturer’s claims for the Brush Recorder are that hysteresis is 
limited to + 1 per cent of the full-scale value. Individual channel 
modules were calibrated with their respective sensors and associated 
equipment before final placement. Calibration figures were prepared 
at five volts excitation using a plus and minus five-foot dynamic head 
pressure fluctuation achieved by lowering and raising the sensors in a 
sea-water tank. After numerous trials, a linear variation of + 1-1/2 
per cent from best straight line was realized for each transducer. 

The harbor data obtained using the Snodgrass Mark IX sensor was 
recorded on an Esterline-Angus model AW strip-chart recorder. Attached 
power supply and bridge circuit components were part of the general 
purpose Mark IX shore wave recorder system. Wave data were recorded 
from sensor No. 3 (harbor) with a curvilinear trace presentation. The 
chart speed utilized was 3/4 inch per minute (0.29 millimeter per second). 

Calibration circuitry unique to the Mark IX shore wave recorder is 
incorporated as a functional component of the system. Calibration 
checks were routinely taken and adjustments were made for deviations 
from set current values for full-scale deflection. 
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V. ANALYSIS OF RESULTS 



A. INTERPRETATION OF POWER SPECTRA 
1. Measured Power Spectra 

A power spectrum of a stationary random process is considered to 
be the statistical estimate of the energy according to frequency bands 
obtained through integration of the spectral density. This spectral 
representation is the Fourier transform of the auto-covariance function. 
This concept of spectral analysis is based on an assumption that the 
energy contained in wind waves and swell can be represented through 
linear theory by the sum of the energy contained in many individual wave 
trains. Therefore, the total area under the spectral density curve will 
be proportional to the total energy. 

The power spectra traces for the three sensors for six sets of 
wave measurements are plotted in Figures 12, 16, 20, 24, 28 and 32. The 
individual traces for all three sensors are superimposed over the same 
frequency interval scale, but the scale of the spectral density ordinate 
for the seaward sensors (No. 1. and No. 2) and the harbor sensor (No. 3) 
is ten to one, respectively. The power spectrum for each sensor was 
computed from time records digitized at 0.63 second intervals giving a 
Nyquist frequency of 1.26 seconds, a resolution of 0.00389 Hz, and a 
total effective record length of 21.5 minutes. Due to the increasingly 
significant hydrodynamic attenuation factors with increasing frequency 
(e.g., 0.40 at 0.15 Hz), the high frequency range of the spectral plots 
was limited to 0.15 Hz (6.7 sec). 

Sensors No. 1, No. 2, and No. 3 were hydraulically high-pass 
filtered at time constants of 47, 50, and 60 seconds, respectively. This 
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FIGURE 12 



Measured Power Spectra, High Tide, 17 Dec 70 (1330) 




Frequency (Hertz) 



FIGURE 13 



Phase Spectrum, High Tide, 17 Dec 70 (1330) 
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FIGURE 14 



Incident, Reflected, and Transmitted Power Spectra 
High Tide, 17 December 1970 (1330) 




Coherence Spectrum, High Tide, 17 December 1970 (1330) 
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FIGURE 16 



Measured Power Spectra, High Tide, 1 Dec 70 (0800) 




FIGURE 17 



Phase Spectrum, High Tide, 1 Dec 70 (0800) 
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FIGURE 18 



Incident, Reflected, and Transmitted Power Spectra 
High Tide, 1 December 1970 (0800) 




Coherence Spectrum, High Tide, 1 December 1970 (0800) 
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FIGURE 20 



Measured Power Spectra, Low Tide, 1 Dec 70 (2000) 
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EIGURE 21 

Phase Spectrum, Low Tide, 1 Dec 70 (2000) 
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FIGURE 22 

Incident, Reflected, and Transmitted Power Spectra 
Low Tide, 1 December 1970 (2000) 




FIGURE 23 

Coherence Spectrum, Low Tide, 1 December 1970 (2000) 
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FIGURE 24 



Measured Power Spectra, Low Tide, 1 Dec 70 (2330) 
<u 
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FIGURE 25 

Phase Spectrum, Low Tide, 1 Dec 70 (2330) 
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FIGURE 26 



Incident, Reflected, and Transmitted Power Spectra 
Low Tide, 1 December 1970 (2330) 




Coherence Spectrum, Low Tide, 1 December 1970 (2330) 
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FIGURE 28 



Measured Power Spectra, Low Tide, 2 Dec 70 (2030) 
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FIGURE 29 

Phase Spectrum, Low Tide, 2 Dec 70 (2030) 
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FIGURE 30 



Incident, Reflected, and Transmitted Power Spectra 
Low Tide, 2 December 1970 (2030) 




Coherence Spectrum, Low Tide, 2 December 1970 (20 30) 
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FIGURE 32 



Measured Power Spectra, Low Tide, 14 Dec 70 (1600) 
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Phase Spectrum, Low Tide, 14 Dec 70 (1600) 
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FIGURE 34 



Incident, Reflected, and Transmitted Power Spectra 
Low Tide, 14 December 1970 (1600) 




Coherence Spectrum, Low Tide, 14 December 1970 (1600) 
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caused substantial dampening (filtering) at an exponentially increasing 
rate for spectral components below the frequencies corresponding to the 
respective time constants. Taking the dampening into consideration, 
tangible spectral energy density was still evident at the lower frequencies 
dealt with below 0.02 Hz (50 sec). These low frequency, filtered, spectral 
densities have been included in the power spectrum plots for information 
with regard to their relative, rather than absolute, spectral density 
values. Of particular note is the strikingly similar energy density 
spectra of seaward sensors No. 1 and No. 2 (both similarly damped) in 
this long wave, low frequency range. Considering the relatively short 
distance between the two sensors and the long wave lengths for the spectral 
wave components within this range of nearly equal spectral densities, it 
is apparent that the instrument calibrations are consistent and accurate. 

All of the spectra analyzed reveal prominant peaks within the 
period range from 50 to 60 seconds; sensor No. 3 (harbor) in the vicinity 
of 50 seconds and sensors No. 1 and No. 2 (seaward) nearer 60 seconds. 
Confirmation of the spectral prominances for these sensors in the 50 to 
60 seconds range was also found in earlier studies. In particular, 

Hudson [1949] indicates persistent spectral peaks at 55 seconds in 
Monterey Harbor. Wilson, Hendrikson, and Kilmer [1965] found through 
the technique of residuation analysis of long-wave records that spectral 
peaks occur frequently at both 60 and 45 seconds for waves within and 
incident at the mcuth of the harbor. 

There is a remarkably high, narrow-banded spectral peak for 
sensor No. 1 centered at a frequency corresponding approximately to a 
14-second period in all the wave records analyzed. This persistent 
feature is particularly interesting in that it coincides with the 
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theoretical anti-nodal line prediction of maximum wave amplitude for the 
prearranged sensor location (a 14-second period quarter wave length 
separation between sensor No. 1 and No. 2, and sensor No. 2 and the 
breakwater face at MLLW) . Correspondingly, a relatively small and flat 
energy density spectrum at sensor No. 2 is found which is located at the 
theoretical nodal position for a 14-second period wave. Spectrum 
analysis thus provides confirmation that a fairly well-defined partial 
standing wave phenomenon exists to seaward of the reflecting breakwater. 
Standing waves were at times observed distinctly in the wave records 
when a well defined 14-second period sinusoidal swell was recorded at 
sensor No. 1, while at the same time a relatively flat nodal area was 
observed at sensor No. 2 (Figure 36). 

In all six wave recordings, sensor No. 1 shows appreciable 
energy in the special frequency range of interest, 0.11 to 0.044 Hz 
(9 to 23 sec), whereas, sensor No. 2 shows considerably less energy but 
does have parallel persistency in spectral peaks for approximately the 
same frequency range. The subsidiary prominant spectral peaks that appear 
at frequencies above 0.11 Hz (9 sec) will be interpreted in the following 
section. 

The power spectra for sensor No. 3 (harbor) is for the most part 
comparatively flat in the frequency range of special interest (0.11 to 
0.044 Hz). In contrast to the high peaks and deep valleys of the seaward 
spectra, the harbor spectrum drops off to significantly low values above 
0.11 Hz (9 sec) and has generally persistent high values below 0.044 Hz 
(23 sec). The spectral peak in the low frequency range closely parallels 
the seaward spectra in both persistence and frequency. This increase in 
spectral amplitudes in the lower frequencies is partially obscured due to 
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Analog Trace Illustrating Partial Standing Wave* Phenomenon 
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the inherent hydraulic filtering (60-second time constant) . Although 
this feature does not prove a serious hindrance for the purpose of this 
study, the filtering does limit quantitative interpretation of the trans^ 
mitted low frequency energy. The definite increase in energy transmission 
toward the lower frequencies can partly be explained through the generali- 
zation that the breakwater acts as a moderately selective filter of short 
wave energy transmission. 

2. Incident and Reflected Sea Surface Power Spectra 

The incident and reflected sea surface power spectra are 
theoretically derived quantities already discussed in this paper under 
theoretical considerations. The incident and reflected spectral wave 
amplitudes are computed for respective spectral frequency components. 

From these analytically developed spectral amplitudes the power spectra 
were computed and are plotted in Figures 14, 18, 22, 26, 30 and 34. 

The sea and swell wave energy is seen to be concentrated in the 
frequency range, 0.11 to 0.044 Hz (9 to 23 sec), for the incident and 
reflected wave components, as expected. The transmitted wave spectrum 
is the same as the No. 3 spectrum, and is also plotted in the above 
figures . 

The one prominant spectral feature apparent in both the incident 
and reflected wave spectra which is centered at the 0.128 Hz (7.8 sec) 
frequency band is an artificial energy spike. This anomaly can be 
readily explained as the result of nodal lines occurring at both seaward 
sensors (No. 1 and No. 2) at approximately the same wave period. The 
two seaward sensors are so positioned that they theoretically measure 
minimum (nodal) pressure fluctuations in the vicinity of this spectral 
frequency band (roughly 7 to 8 seconds) . The separation of the two wave 
sensors corresponds to a half-wave length of a 7.5 second wave period. 
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This unavoidable consequence caused by the partial standing wave 
creates the lack of discernible pressure fluctuation measurements bet- 
ween the two seaward sensors for this double nodal situation. The 
spectral densities in the vicinity of this frequency band are theoretic- 
ally indeterminate, in that the measurements give an extremely high 
peaked, meaningless energy spike. This undesirable spectral feature 
limits the high frequency range of interest for this study to 0.11 Hz 
(9 sec) . It would have been necessary to either position the sensors 
differently (larger interval spacing) or use three sensors seaward of 
the breakwater to have resolved these frequency bands. 

B. INTERPRETATION OF CROSS SPECTRA 

The cross spectrum is the Fourier transform of the cross-covariance 
function; the cross-covariance between two time-series records describes 
the general dependence of the values of one time-series record on the 
other. The cross-covariance function has essentially the same form as 
the auto-covariance function except that one independent time-series 
record, X^, is examined (lagged) through another time-series record, X^, 
rather than through itself as with the case of the auto-covariance 
function. The mutual dependence between the two time-series records can 
be quantitatively expressed through the cross spectra phase and the 
coherence function. The latter, being an estimator of the relative phase 
stability between the two cross-correlated records, is represented by a 
correlation coefficient from zero to one. When the coherence is near zero, 
the two records are said to be uncorrelated, or incoherent. The correlation, 
or degree of coupling between two random records, increases as the coherence 
approaches unity. 
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Coherence spectra curves for each set of wave observations are 
shown in Figures 15, 19, 23, 27, 31, and 35. Low coherence or non- 
correlation between any two random records of the phenomenon being 
studied can be due to any number of causes of which the more creditable are 
as follows : 

1. Incident wave crests not parallel to the reflecting surface, 
thus presenting a complex, multi-peaked, incident-reflected wave surface. 

2. Significant scattering of the incident wave energy vice normal 
reflection due to the rough angular reflecting surface of the rubble-mound 
breakwater (predominately high frequency reflection) . 

3. Coherent energy levels so low that they are lost in background 

noise. 

Phase spectra curves for each set of wave observations are shown in 
Figures 13, 17, 21, 25, 29, and 33. Based on standing wave theory, for 
waves normally incident, the phase relationships between seaward sensors 
No. 1 and No. 2 should be near zero degrees for the low frequencies 
approaching tidal period; near 90 degrees for frequencies of waves where 
the separation between the seaward sensors approaches one-quarter wave 
length; and near 180 degrees for one-half wave length separation. This 
theoretical curve is described by the dashed curve superimposed on the 
measured phase difference curve in each figure. The comparison between 
theory and data is somewhat complex and will be interpreted for different 
ranges of the frequency spectrum. 

The frequency range from zero to 0.044 Hz (23 sec) has a theoretical 
phase difference varying from zero to 25 degrees, increasing with 
increasing frequency. The phase curves compare reasonably well on all 
the cross spectra phase plots in this frequency range. The coherence is 
high (average 0.90). 
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In the frequency range 0.044 to 0.05 Hz (23 to 20 sec) the large 
discrepancy in the phase curves for the theoretical and measured phase 
spectra may be explained by the fact that it is a spectral region of 
poor coherence (average 0.20). A plausible reason for this lack of 
correlation is simply the absence of any energy in the frequency band 
(note the minor energy density in this frequency range) . 

The phase spectra in the frequency range from 0.05 to 0.15 Hz (20 to 
6.7 sec) shows variable deviations between the theoretical and the 
measured curves. The cross-correlated phase data associated with high 
tide on 1 December and 17 December 1970 (Figures 13 and 17) compare 
within a maximum deviation of approximately 25 degrees. The other cross- 
correlated records show somewhat less correspondence. All of the records 
tend toward finer correspondence of the measured and theoretical phase 
differences with increasing frequency. Coherence drops almost linearly 
from moderate (average 0.70) to poor (average 0.20) in this range. For 
two wave records which tend toward negligible energy density in the 
vicinity of 0.10 Hz (10 sec) (Figures 12 and 24), there is an expected 
although slight dip in coherence. Reasons for the decreasing coherence 
tendency within this frequency range are uncertain. The most plausible 
reason is the tendency for increased scattering as wavelength decreases 
with increasing frequency. The reflecting surface area of the angular 
shaped armor units relative to the decreasing wave length assumes a 
controlling factor which probably magnifies the scattering effect 
(oblique incidence to the reflecting surface) . Another possible 
explanation is oblique wave incidence of multidirectional free wave 
components of the same frequency but of different origin. 
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C. TRANSMISSION, REFLECTION, AND ENERGY DISSIPATION COEFFICIENTS 

Transmission, reflection, and energy dissipation coefficients are 
dimensionless parameters which describe the usefulness of a permeable 
breakwater. Transmission and reflection coefficients computed from the 
six selected time-series records (Table I) are presented according to 
respective spectral frequency bands in Figures 37 and 41, and the energy 
dissipation coefficients are presented in Figures 38 and 42. For all 
measurements, winds were nominal and the direction of sea and swell was 
visually observed to be within ten degrees of normal incidence to the 
breakwater. Data scatter is significant only at the high and low 
boundaries of the spectral frequency range of interest, 0.11 and 0.044 Hz 
(9 and 23 sec) . The reasons for the scatter are directly related to the 
previous discussion of power and cross spectra pertaining to the anomalous 
measurements in these select border frequency bands. 

Figures 37 and 38 describe the coefficients at high tide conditions 
(average four feet MLLW) for two time-series records (Table I). Both 
wave records were at incipient wave overtopping of the breakwater 
(Figures 39 and 40) . The coefficient of transmission has a nearly 
uniform value of about 0.10 from 0.11 to 0.067 Hz (9 to 15 sec), below 
which increases almost exponentially as frequency decreases. This 
increase can be aptly explained as the result of the longer wave length 
of the incident lower frequency wave components. 

The wave reflection coefficient decreases almost linearly with 
increasing frequency due to the progressively dominating influences of 
the scattering of the incident wave energy as the wave length decreases. 
The energy dissipation coefficient, which is primarily due to the random 
scattering, turbulence, and the internal frictional resistance of the 
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Monterey Breakwater Wave Data from Field Study 
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Note: (1) Condition of incipient wave 

high tide conditions. 
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Coefficients of Reflection and Transmission 



FIGURE 37 

Coefficients of Reflection and Transmission, High Tide Conditions 
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Frequency (Hertz) 




Coefficients of Energy Dissipation 



FIGURE 38 

Coefficients of Energy Dissipation, High. Tide Conditions 
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Frequency (Hertz) 




FIGURE 39 



Breakwater at High Tide, 1 December 1970 
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FIGURE 40 



Breakwater at High Tide, 17 December 1970 
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rubble-mound structure, Increases consistently with an increase in 
frequency, noticeably counter to the decreasing reflection coefficient. 

The reflection coefficient appears clearly related to the frequency 
variable. The coefficient of reflection occurring at the frequency of 
peak energy concentration, 0.071 Hz (14 sec) is approximately 0.50. The 
shape of the best-fit curve to the transmission coefficient values does 
not lend itself to a meaningful average value, although a coefficient of 
0*10 at the frequency of peak energy concentration is representative of 
the sea and swell energy. A value of approximately 0.70 occurs for the 
energy dissipation coefficient at the frequency of peak energy 
concentration . 

Figures 41 and 42 describe the coefficients at low tide conditions 
(average one foot MLLW) for four selected time-series records (Table I) . 
No wave overtopping of the breakwater occurred on any of these occasions 
(Figure 43) . Interpretation of these records is not as clear cut as 
for the high tide conditions. Data scatter was somewhat more prevalent 
and not as closely confined to the boundaries of the spectral frequency 
range of interest, 0.11 Hz (9 sec) and 0.044 Hz (23 sec). 

The coefficient of reflection may be seen to reach a maxima at a 
frequency of 0.08 Hz. A reason as to this non-linear reflection curve 
(Figure 41) is a direct consequence of the structural condition of the 
breakwater. The exposed seaward slope has altered its configuration to 
a small degree and not uniformly. Rubble structures over the years will 
follow a process of disintegration; that is, wearing away, or dislodging 
stone by stone, rather than substantial collapse. This pseudo altered 
structure actually forms a more stable base. The portion of the break- 
water more frequently exposed to attack by high breaking waves and the 
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FIGURE 41 

Coefficients of Reflection and Transmission, Low Tide Conditions 
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FIGURE 42 

Coefficients of Energy Dissipation, Low Tide Conditions 
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Frequency (Hertz) 





FIGURE 43 

Breakwater at Low Tide, 14 December 1970 
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surge through the permeable structure lies in the upper section of the 
exposed breakwater. It is evident that the Monterey Breakwater has 
had armor stones dislodge from the upper section and subsequently 
settle along the breakwater slope in the lower section of the exposed 
breakwater. This has resulted in a decrease in slope on the seaward 
side to provide less verticality to incident wave energy. At low tide 
conditions, the higher energy wave components near the 14-second period 
band are incident to a greater extent on the more vertical section of 
the breakwater face higher up the slope and subsequently have a higher 
reflection coefficient. The lower energy wave components on both sides 
of the period of maximum energy concentration are dissipated and scattered 
to a larger extent by the greater number of less vertical angular 
reflecting surfaces along the direction of sea and swell incidence. The 
higher frequency energy components are scattered to an even greater 
degree due to their shorter wave lengths with considerable reduction in 
reflection coefficient. 

A reflection coefficient of approximately 0.60 occurs at the 
frequency of peak energy concentration, 0.071 Hz (14 sec), and is 
slightly higher than the high tide conditions. This is presumably due 
to the absence of any wave overtopping as well as to the unique reflective 
considerations just discussed. The coefficient of transmission is 
similar in trend but a little lower in magnitude compared to the high 
tide conditions, having a constant coefficient of approximately 0.03 
at frequencies above about 0.06 Hz. The smaller value of the coefficient 
is due to the denser packing of the rubble units in the lower section of 
the breakwater. The energy dissipation coefficient is similar to that 
for the high tide conditions in that it is also inversely related to 
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the reflection coefficient* A value of approximately 0.60 for the 
energy dissipation coefficient occurs at the frequency of peak concen- 
tration of energy, 0.071 Hz (14 sec), and is slightly lower than the 
high tide conditions. 

D. COMPARISON WITH MODEL STUDY 

Tangible results regarding transmission characteristics through 
model permeable rubble-mound structures in past studies have generally 
been derived for monochromatic waves generated in flume-type wave tanks 
or extensive harbor models. These small-scale models are designed and 
operated, and the test results are converted to prototype dimensions on 
the basis of transference equations derived from Froude's model laws. 

In addition, to preserve accurate results, the model breakwaters are 
constructed geometrically similar to the prototype. Harbor model 
studies are generally built with distorted scales which have adverse 
effects on the accuracy of test results when subjected to sea and swell 
[Hudson and Moore, 1950]; therefore, wave transmission studies are 
carried out primarily in wave tanks or flumes. 

A model test of Monterey Harbor was constructed on a l:40-scale 
breakwater structure at the U. S. Army Engineer Waterways Experiment 
Station to determine wave transmission characteristics of a proposed 
rubble-mound breakwater [Davidson, 1969]. Measured heights of trans- 
mitted and incident waves for a simulated 15-second monochromatic wave 
are tabulated in Table II. The results of the wave transmission test 
show that in general, heights of the transmitted waves were from 25 to 
40 per cent of those of the incident waves for waves up to incipient 
overtopping of the structure. For purposes of comparison, the model 
breakwater differed from the prototype rubble-mound structure in that 
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TABLE II 



Monterey Harbor Wave Transmission Data 
from Hydraulic Model Investigation 
(after Davidson, 1969) 



H I 

INCIDENT 
WAVE HEIGHT 
(ft) 


TRANSMITTED WAVE HEIGHT (ft) 
MEASURED AT: 


L/4 


L/2 


3.0 


0.8 


0.8 


4.5 


1.2 


1.2 


6.0 


1.6 


1.2 


7.5 


2.0 


1.6 


9.0 


2.4 


2.8 


10.5 


4.0 


4.0 



Note: (1) Hj. is the height of the incident wave 

measured at the structure site before 
the structure was installed in the model. 



(2) Significant overtopping occurred with 
H^ = 10.5 feet. 

(3) Monochromatic wave period was 15 seconds (simulated). 

(4) H^ was measured at distances of one-quarter 

wave length (L/4) and one-half wave length 
(L/2) . 

(5) Wave flume used was 119 feet long, 5 feet 
wide, and 4 feet deep. 
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it had a crest height of plus 14 feet MLLW (four feet higher) with 
12-ton tribars for the outer cover layer. Transmitted wave heights, 
plotted in Figure 44, were measured at distances of one-half and one- 
quarter of a 15-second wave length from the breakwater centerline. 

The prototype data in Table I are plotted in Figure 44 to offer a 
comparison between the field and hydraulic model measurements. The 
prototype breakwater has a crest height of ten feet MLLW, which accounts 
for near wave overtopping conditions at lower incident wave heights 
than for the hydraulic model study. The heights of transmitted waves 
were from 7 to 12 per cent of the incident waves for conditions up to 
and including incipient wave overtopping of the structure. 

Diffraction effects for the prototype breakwater have not been 
taken into account up to this point in the analysis. During the 
periods of field measurement no waves due to diffraction were visually 
observed in the vicinity of sensor No. 3 (harbor). Appropriate dif- 
fraction diagrams designed for semi-infinite impermeable breakwaters 
indicate a diffraction coefficient of less than 0.10 for a 15-second 
wave at the location of sensor No. 3 [Corps of Engineers, 1961]. 
Laboratory diffraction tests have shown that the theoretical values 
of diffracted wave heights within the geometric shadow of a breakwater 
exceed the actual values by about ten per cent [Dunham, 1950]. It is 
the author’s opinion, in consideration of the fact that the geographic 
location of the sensor is far within the geometric shadow and at a 
sharp angle from the breakwater tip (500 feet from the tip of the break- 
water and approximately 80 degrees from the incident wave direction) , 
that the validity of diffraction theory is suspect at the location of 
the harbor sensor. The diffraction coefficient is considered to be well 
under 0.10, probably approaching zero over the period range dealt with. 
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□ MODEL TEST DATA 
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FIGURE 44 

Transmitted Wave Heights for Protopype Study 
and Hydraulic Model Investigation 
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INCIDENT WAVE HEIGHT «») 



As a consequence of discounting diffraction as contributing to 
the wave heights measured at the harbor sensor, the coefficients of 
transmission from the model studies appear even more conservative, 

E. ERROR ANALYSIS 

The errors encountered in this field study can be classified as 
either data-reduction, instrument location, or numerical-analysis 
oriented. Each of these error possibilities are examined below in 
order to determine their effect on the calculated results of this study. 

1. Data Reduction 

Time-induced errors may arise from synchronization of the wave 
records which is critical to spectral phase computation. The recording 
of the raw wave data from the two seaward sensors on analog traces was 
simultaneously done on a multi-channel recording instrument , which 
insures exact synchronization. 

The harbor sensor (No. 3) was recorded as a curvilinear trace 
but digitized by a Calma 480 digitizer which accepts only rectilinear 
coordinates. The digitized data were subsequently reconverted to 
curvilinear coordinates. Linear interpolation was then used to match 
data points with the two seaward (sensors No. 1 and No. 2) records. 

The sampling interval was reduced from 0.80 to 0.63 seconds which 
maintained high resolution. Reduction and conversion of data in this 
manner resulted in very minor raw data errors . 

A relation, obtained from linear wave theory, between wave 
height and recorded differential pressure was used to adjust the 
observed wave height according to wave frequency and sensor depth. 

Recent research by Bergan, Torum, and Traetteberg [1968] involving 
a comparison of wave measurements by a pressure-type wave gauge and 
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a continuous wire wave gauge found that the height difference of 
individual waves observed is significant. However, this same study 
confirms that the analysis of pressure gauge records of random sea 
waves gives reasonably accurate wave spectra and wave-height distribu- 
tions. The agreement between significant wave heights obtained from 
the pressure and continuous wire gauge was within ten per cent. For 
the water depths and wave periods dealt with in this study, the error 
is considered well within ten per cent. 

2. Instrument Location 

The accuracy achieved in sensor positioning by transit and 
range line is considered within two feet. A bottom survey in the 
vicinity of the breakwater (Figure 4) was conducted by the author 
utilizing a portable fathometer with positioning fixed by transit, 
measuring tapes, and range lines. Depth accuracy was established 
within one foot. Therefore, errors in sensor positioning and bottom 
depth are considered to have negligible effect on the calculated 
results . 

3. Numerical Analysis 

The effect of approximating a finite time-series record with 
a step function and transforming the step function through computer 
spectral analysis can result in a number of misleading interpretations 
if not handled properly. Initially, this method of approximation by a 
step function can have considerable effect due to aliasing. Due to the 
high frequency cut-off and small sampling intervals of this study, 
aliasing, or the superposition of high frequency energy on the recorded 
wave form, is negligible. 

Another source of numerical error is caused by the effect of 
transforming the record into discrete values of energy-density as a 
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function of discrete frequency bands. The energy is forced into 
certain band widths and centered about the middle frequency. In an 
effort to develop independent spectral amplitudes and modify any 
errors introduced, the records were smoothed by applying a Parzen 
tapered lag window to the auto-co variance function. 

The assumption of an infinitely long record in theory and the 
use of a finite record in the analysis produces an error which can be 
expressed by confidence limits of the computed spectrum. The confidence 
limit is the probability that the true spectrum lies between the actual 
spectrum estimated from a finite sample. These bounds are a function 
of the number of degrees of freedom of the wave record, which is a 
function of record length and number of lags. To achieve the necessary 
balance between adequate resolution and reliable spectral density esti- 
mates, the number of lags was established at ten per cent of the total 
number of data points digitized (205 lags). This limited the number of 
degrees of freedom in each spectrum to 20 for each 21.5-minute time 
series. This, in turn, determined the confidence limits at the esti- 
mated power spectrum. For example, the 80 per cent confidence limits 
for 20 degrees of freedom are interpreted as meaning that one can state 
with 80 per cent accuracy that the true spectral estimate is greater 
than 0.71 but less than 1.60 at the actual power spectrum, as the 
spectrum is distributed according to chi-square distribution [Jenkins 
and Watts, 1969]. 
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VI. CONCLUSIONS AND SUMMARY 



The conclusions drawn from this field study are summarized as 
follows : 

A. The method of resolving incident and reflected wave heights 
for a partial standing wave through the use of two fixed wave sensors 
was for the most part successful. The basic drawback was the inability 
to perform high frequency analysis in the period range of nine seconds 
and less. If the sensors had been positioned with a greater separation, 
the double nodal condition would have been displaced toward higher 
frequencies. The use of three sensors instead of two to seaward of 

the breakwater w T ould have served the same purpose to resolve the 
incident and reflected waves in frequency bands above nine seconds. 

B. The treatment of the wave transmission problem as a linear, 
stationary, random process is considered a reliable and useful tech- 
nique. Power spectra complimented with cross spectral correlation 
resolved the wave field to seaward of the breakwater into incident and 
reflected spectral components. Confidence bounds of 0.71 and 1.60 for 
a confidence interval of 80 per cent were realized using a record 
length of 2048 data points with a maximum number of 205 lags. 

Increased record length with a decrease in the ratio of number of 
lags to data points would have satisfied the requirements of good 
resolution and improved the confidence bounds. Whereas the sampling 
interval was fixed, the record length and number of lags were 
arbitrarily determined. 

Coherence was lower in the high frequency domain, and is considered 
a result of the multi-directional scattering effect of incident waves 
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by the rough, angular, reflecting surfaces of the rubble-mound 
breakwater, 

C. The transmission, reflection, and energy dissipation 
coefficients calculated from the prototype data followed the trends 
generally expected for wave energy incident on a rubble-mound break- 
water. For high tide conditions, the reflection coefficient decreased 
linearly with increasing frequency. The energy-dissipation coefficient 
increased with increasing frequency in contrast to the reflection 
coefficient. The transmission coefficient was found to be approxi- 
mately constant from the high frequency cut-off of 0.11 Hz (9 sec) to 
about 0.067 Hz (15 sec) where it increased abruptly with decreasing 
frequency. Coefficients of 0.10 for transmission, 0.50 for reflection, 
and 0.70 for energy dissipation occurred at the frequency of peak 
energy concentration, 0.071 Hz (14 sec). At low tide conditions, the 
coefficients were slightly different. The approximate reflection 
coefficient was slightly higher; the transmission coefficient lower, 
and the energy dissipation coefficients less than the high tide 
conditions . 

D. The numerical results obtained from the prototype transmission 
tests and from the proposed Monterey Breakwater model did not agree 
particularly well. For the hydraulic model study, the heights of 
transmitted waves were reported to be from 25 to 40 per cent of those 
of the incident waves, whereas the prototype field measurements were 
from 7 to 12 per cent within the same general frequency range and with 
no significant wave overtopping. It is important to note that signifi- 
cant wave height is the height parameter used in the prototype trans- 
mission study; whereas, the actual heights of the monochromatic waves 
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were used in the hydraulic laboratory investigation. The possibility 
remains that the model laws of similarity for prototype scaling and 
the effect of artificially generated waves produce over-estimated 
results. It may be concluded from this study of the Monterey Harbor 
breakwater that the design of rubble-mound breakwaters in general, 
based on model studies of wave energy transmission, may be conservative. 
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structure in the natural environment. Transmission has been studied 
heretofore exclusively with scale models. The coefficients of trans- 
mission are considerably less for the prototype than predicted by a 
Corps of Engineers model study of the Monterey Breakwater. The differ- 
ences are apparently related to more wave energy scattering and dissi- 
pation due to turbulence than predicted by the model. 
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